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Abstract 

Severe acute respiratory syndrome (SARS) is a respiratory disease caused by a newly found virus, called SARS coronavirus. In 
this study, the cleavage mechanism of the SARS coronavirus main proteinase (M pro or 3CL pro ) on the octapeptide NH 2 - 
AVLQ i SGFR-COOH was investigated using molecular mechanics and quantum mechanics simulations based on the experimental 
structure of the proteinase. It has been observed that the catalytic dyad (His-41/Cys-145) site between domains I and II attracts the n 
electron density from the peptide bond Gln-Ser, increasing the positive charge on C(CO) of Gin and the negative charge on N(NH) 
of Ser, so as to weaken the Gln-Ser peptide bond. The catalytic functional group is the imidazole group of His-41 and the S in Cys- 
145. N S1 on the imidazole ring plays the acid-base catalytic role. Based on the “distorted key theory” [K.C. Chou, Anal. Biochem. 
233 (1996) 1-14], the possibility to convert the octapeptide to a competent inhibitor has been studied. It has been found that the 
chemical bond between Gin and Ser will become much stronger and no longer cleavable by the SARS enzyme after either changing 
the carbonyl group CO of Gin to CH 2 or CF 2 or changing the NH of Ser to CH 2 or CF 2 . The octapeptide thus modified might 
become an effective inhibitor or a potential drug candidate against SARS. 

© 2004 Elsevier Inc. All rights reserved. 
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A new coronavirus (SARS CoV) * 1 has been found as 
the etiological agent of the disease SARS [1-4]. SARS 
CoV is a positive-sense, single-stranded RNA virus with 
the largest viral RNA genome known to date [5,6]. The 
SARS CoV genome is composed of around 29,700 
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nucleotides. The replicase gene alone encompasses more 
than 21,000 nucleotides, 2/3 of the whole genome se¬ 
quence. The replicase gene encodes two overlapping 
polyproteins, ppla (4377 a.a.) and pplab (7072 a.a.). 
SARS CoV main proteinase (M pro ) cleaves the polypro¬ 
tein at no less than 11 conserved sites, a process initiated 
by the enzyme’s own autolytic cleavage from ppla and 
pplab [7,8]. The released polypeptides mediate all of 
the functions required for SARS viral replication and 
transcription. The functional importance of the M pro 
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in the viral life cycle makes it an attractive target for the 
drugs design against SARS [7-10]. 

Anand et al. [7] suggested that the rhino virus 3C pro 
inhibitor AG7088 could well serve as a starting point 
for an anti-SARS drug based on the theoretical homol¬ 
ogy model of SARS CoV M pro . Chou et al. [9] found the 
fitting problem of AG7088 to the binding pocket of 
SARS CoV M pro , and they suggested its derivative 
KZ7088 as a better starting point. Sirois et al. [11] did 
virtual screening for SARS-CoV protease based on 
KZ7088 pharmacophore points. Chou et al. [9] also sug¬ 
gested an octapeptide NH 2 -AVLQ 1 SGFR-COOH that 
can be converted to an inhibitor of SARS CoV M pro by 
some chemical modification. This octapeptide was syn¬ 
thesized, proved to be a cleavable peptide, and showed 
strong inhibiting activity in a tube test [12]. Du et al. 
[13] analyzed the cleavage mechanism of the octapeptide 
by SARS CoV M pro and discussed the possible chemical 
modification of the octapeptide using molecular dynam¬ 
ical and quantum mechanical simulations based on the 
homology structure of SARS CoV M pro , to convert it 
into a stable and effective drug. This progress on finding 
inhibitors against the SARS enzyme is also summarized 
in a recent review [10]. 

The experimental structure of SARS CoV M pro was 
first determined by Rao’s research group [8]. They pro¬ 
vided several crystal structures of SARS CoV M pro at 
different pH conditions and the complex structure of 
SARS CoV M pro with ligand hexapeptidyl CMK [7,8]. 
The experimental structure and the theoretical homol¬ 
ogy structure of SARS CoV M pro share many structural 
similarities, such as the catalytic active region and the 
Cys-His cleavage dyad [7,8]. However, there are several 
differences between the theoretical and the experimental 
structures. For examples, the experimental structure of 
SARS CoV M pro is a dimer, the protomer A has cata¬ 
lytic activity, and the protomer B has no activity. How¬ 
ever, the N finger of protomer B, which plugs in the 
active cleft of protomer A, plays an important role in 
both dimerization and maintenance of the active form 
of the enzyme [8]. These new findings from the experi¬ 
mental structure of SARS CoV M pro provide a reliable 
basis for inhibitor design. 

Method 

The octapeptide NH 2 -AVLQ I SGFR-COOH shows 
strong inhibiting activity in vitro [12]. However, the use 
of native peptide for clinical applications has been ham¬ 
pered by their intrinsic metabolic instability (indigestion 
by enzymes) and poor absorption through membrane. A 
cleavable peptide of SARS CoV M pro means that it pos¬ 
sesses a competent binding with its receptor and a 
breakable peptide bond in the cleavage site. To convert 
a cleavable peptide into a stable and effective inhibitor, 


the peptide has to be modified according to the “distort¬ 
ed key” theory [14,15]. One important step for the chem¬ 
ical modification is a good understanding of the 
cleavage mechanism, particularly the changes in chemi¬ 
cal bonds [16-22]. 

The protease-susceptible sites in a given protein or 
peptide usually extend to an octapeptide region [14,15]. 
The corresponding amino acid residues are sequentially 
symbolized by eight subsites R 4 , R 3 , R 2 , Ri, Ry, R 2 s 
R 3 /, R 4 / and the eight matching positions in the protease 
by S 4 , S 3 , S 2 , Si, Si', S 2 /, S 3 /, S 4 /, respectively (see Fig. 3 
of [15]). Occasionally, the susceptible sites in some pro¬ 
teins may contain one subsite less or one more. How¬ 
ever, eight amino acid residues are the most common 
case. Although the protein being cleaved contains much 
more than eight amino acid residues, only an octapep¬ 
tide segment fits in the active region of the protease. 
Therefore, our research will focus on the cleavability 
of octapeptides. 

As shown in Fig. 1, the cleavage point is always on 
the peptide bond between Ri and R^. In Fig. 1A we 
use the combination of a thin line and a dashed line to 
represent the conjugate n property in the peptide bond. 
In Fig. IB the breakable peptide bond between Ri and 
Ri is converted to a strong “hybrid peptide bond” 
through a chemical modification and hence cannot be 
cleaved by the protease. The modified octapeptide can 
still bind to the active site, but it cannot be cleaved by 
the enzyme, just like a distorted key that can neither 
open the lock nor be pulled out. The octapeptide thus 
modified will automatically become a desired inhibitor 
candidate. 

The octapeptide AVLQSGFR is the first suggested [9] 
based on the molecular structure of SARS CoV M pro 
and has been proved cleavable by experiments. In this 
research we study the cleavage mechanism, the proper¬ 
ties of the relevant chemical bonds, and the catalytic 
interactions between the octapeptides and the SARS 
CoV M pro using molecular mechanical and quantum 
chemical simulations to provide useful insights for the 
chemical modification. The study was carried out in 
the following four steps: (1) docking the octapeptide to 
SARS CoV M pi ° followed by an energy minimization 
for the complex structure using molecular mechanics, 
(2) computing the atomic charge distribution in the 
binding pocket for the energy-minimized structure using 
ab initio quantum mechanics, (3) computing the molec¬ 
ular energy, chemical bond properties, and atomic 
charges of the octapeptide with the background charge 
distribution [23] of SARS CoV M pro using ab initio 
quantum chemistry, and (4) with the same background 
charge distribution, computing the molecular energy, 
chemical bond properties, and atomic charges of the 
modified octapeptide. The computations have been car¬ 
ried out with the molecular dynamics and standard ab 
initio quantum chemistry at the HF/6-31G* level [23]. 
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Fig. 1. Schematic illustration to show (A) a cleavable octapeptide is chemically effectively bound to the active site of SARS enzyme and (B) although 
still bound to the active site, the peptide has lost its cleavability after its scissile bond was modified from a hybrid peptide bond to a strong bond. 
Adapted from Chou [15] with permission. 


Results 

Docking and energy minimization 

The experimental structure (lujl.pdb) of SARS CoV 
M pro [8] is used in the calculation and shown in Fig. 2A. 
Each protomer structure of the SARS CoV M pi ° dimer 
has three domains. The two protomers A and B form a 
right angle and the join part is in the domain III. The N 


finger of protomer B plugs into the domain II of pro¬ 
tomer A and plays an important role in the catalytic 
activity. The PDB data of SARS CoV M pi ° contains 
300 hydration water molecules, which provide the solva¬ 
tion ion environment with some water molecules joining 
the catalytic reaction. In Fig. 2A the green signs “<” re¬ 
present the solvation water molecules. The docking 
calculation between SARS CoV M pro and the octapep¬ 
tide AVLQSGFR was performed using the molecular 



Fig. 2. (A) Crystal structure of SARS CoV M pro consists of two protomers A and B. Each protomer has three domains and the two protomers form 
a right angle at the domain III. The N finger of protomer B plugs into the domain II of protomer A. The green signs “<” represent the solvation 
water molecules. (B) The complex structure obtained by docking the octapeptide to SARS CoV M pro and a follow-up energy minimization. The 
catalytic dyad of His-41 and Cys-145 is marked by an arrow and the peptide bond to be cleaved of the octapeptide is shown in green. 
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operating environment program package [24]. A total of 
25 docking conformations were obtained, and the one 
with the lowest energy was used for further energy min¬ 
imizing calculation. Fig. 2B shows the complex structure 
of SARS CoV M pro with the ligand octapeptide. The 
positions of catalytic dyad His-41 and Cys-145 are indi¬ 
cated in Fig. 2B. 

According to Anand et al. [7], Yang et al. [8], and 
Chou et al. [9], the catalytic active region is in the pocket 
between domains I and II of protomer A. During calcu¬ 
lation the catalytic active region on which we were 
focusing contains the amino acid residues Cys-22, Gly- 
23, Thr-24, Thr-25, Leu-27, His-41, Val-42, Cys-44, 
Thr-45, Ala-46, Glu-47, Asp-48, Met-49, Leu-50, Asn- 
51, Pro-52, Tyr-54, Phe-140, Leu-141, Asn-142, Gly- 
143, Ser-144, Cys-145, Glu-163, His-164, Met-165, 
Glu-166, His-172, Asp-187, Arg-188, and Gin-189 from 
promoter A and the residues Phe-3, Arg-4, Lys-5, Met- 
6, and Ala-7 from promoter B. But in the practical cal¬ 
culation, the actual region involved contains much more 
than the above-listed amino acid residues. It includes all 

o 

the amino acid residues within 5 A surrounding the 
octapeptide [9,18]. 

Fig. 3A shows the catalytic active region of SARS 
CoV M pro and the location of octapeptide (white). The 
catalytic dyad of His-41 and Cys-145 is shown by a ball 
and stick drawing, the peptide bond Gln-Ser is in green, 
and the N finger in pink. Fig. 3B gives a close look at the 
bridge water molecule (W) between the SARS CoV M pro 
and the octapeptide. The distance 0(W)-C(C0) is 
4.70 A, H^W)- N §1 (His-41) is 2.39 A, H 2 (W)-N(NH) 
is 4.20 A, and H(NH)-S(Cys-145) is 4.12 A. The bridge 


water molecule plays an important role in the proteolytic 
reaction. These results are derived from the complex 
structure obtained by following the aforementioned 
docking and energy-minimization procedures. 

Hydrophobic-hydrophilic interaction 

The hydrophilic and hydrophobic complementary 
interaction is important in the ligand-receptor interac¬ 
tion. According to the a sphere model [25,26], the active 
region of proteinase consists of a series of a spheres. 
Each interaction site between the ligand and receptor 
contains one or more a spheres. There must be at least 
one hydrophobic interaction site. Figs. 4 A and B are 
the hydrophilic and hydrophobic surfaces of octapeptide 
AVLQSGFR and the catalytic active region of SARS 
CoV M pro , respectively. In the figures the green parts 
are hydrophobic surfaces and the blue parts are hydro¬ 
philic surfaces. There is good hydrophilic and hydro- 
phobic complementary interaction between ligand and 
receptor: the hydrophilic parts of the ligand correspond 
to the hydrophilic parts of the receptor and the hydro- 
phobic parts of the ligand correspond to the hydro- 
phobic parts of the receptor. The octapeptide 
NH 2 -AVLQ1SGFR-COOH has four hydrophobic 
amino acid residues: R4(Ala), R3(Val), R3'(Phe), and 
R4'(Arg). Among them R4(Ala) and R3(Val) are ex¬ 
posed in solvent and R3'(Phe) and R4'(Arg) are covered 
by hydrophobic surfaces of the proteinase. The peptide 
bond to be cleaved between Rl(Gln) and Rl'(Ser) is 
just in the join point between the exposed part and the 
covered part. 



Fig. 3. (A) Catalytic active region of SARS CoV M pro and the location of octapeptide (white). The catalytic dyad of His-41 and Cys-145 is shown in 
ball and stick form, the peptide bond Gln-Ser is in green, and the N finger is in pink. (B) Bridge water molecule (W) between SARS CoV M pro and 
octapeptide. The distance Q(W)-C(CO) is 4.70 A, H!(W)-N 51 (His-41) is 2.39 A, H 2 (W)-N(NH) is 4.20 A, and H(NH)-S(Cys-145) is 4.12 A. 
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Fig. 4. (A) Hydrophilic and hydrophobic surfaces of octapeptide NH-AVLQ i SGFR-COOH. The green parts are hydrophobic surfaces and the 
blue parts are hydrophilic surfaces. The four hydrophobic amino acid residues are R4(Ala), R3(Val), R3'(Phe), and R4'(Arg). (B) Hydrophilic and 
hydrophobic surfaces of SARS CoV M pro . The hydrophobic residues R4(Ala) and R3(Val) are exposed in solvent, while R3'(Phe) and R4'(Arg) are 
covered by hydrophobic surfaces of the proteinase. 


Quantum chemical calculations 

The calculations were focused on the effect of SARS 
CoV M pro on the peptide bond to be cleaved. How to 
properly treat a biomacromolecule and the solvent ef¬ 
fects is a nontrivial problem in quantum chemical calcu¬ 
lations. Electrostatic interaction plays the dominant role 
in ligand-receptor combinations and must be taken into 
account in the quantum mechanical study. In view of 
this, we treated the SARS CoV M pro as a background 
with many point charges. The amino acid residues in 
the catalytic region between domains I and II of pro- 
tomer A were divided into seven segments, and their 
atomic charges were separately calculated using the ab 
initio quantum chemical HF/6-31G* method. The 5 
amino acid residues of the N finger of protomer B were 
taken into account also. The 76 amino acid residues in 
eight segments are listed in Table 1. There are 1148 
atomic charges in the background, including all atoms 
in the catalytic cleft; the overlapped atoms are not 
counted. 

Another difficult problem for quantum mechanical 
calculation for a biological system in solution is the sol¬ 
vent effect. In the PDB data [8] of SARS CoV M pro there 


are 300 solvation water molecules. These water mole¬ 
cules not only provide the solvation environment but 
also indicate that some water molecules are involved 
in the catalytic reaction. In our calculation the 300 sol¬ 
vation water molecules were treated as the background 
charges. Each water molecule has 3 point charges. 
Therefore, there are 1148 + 300 x 3 = 2048 background 
point charges. 

In Table 2 we list the Mulliken atomic charges q Mul1 
and electrostatic potential equivalent charges g ESP of 
His-41, which are obtained by fitting atomic charges to 
the electrostatic potential at the van der Waals surface. 
The atomic charges from semiempirical method AMI 
are quite different from the results by ab initio HF/6- 
31G* [23]; in particular; the g ESP charges of AMI do 
not appear reasonable. Because the HF/6-31G* g ESP 
charges reproduce the quantum mechanical electrostatic 
potential on the molecular surface, in this research we 
use q ESP from the HF/6-31G* calculations to illustrate 
our points. The imidazole group of His-41 is a nucleo¬ 
philic agent and its p K = 6.0, near the concentration 
of proton [H + ] in pure water. The imidazole ring is both 
a proton donor and a receptor. It has good catalytic 
activity in the physiological condition [27]. It can be seen 


Table 1 

Distribution of amino acid residues in the active cleft of SARS-CoV M pro : 76 amino acid residues and 1148 atoms are included 


Domain I of A 



Domain II of A 




B 

Al-1 

A 1-2 

A 1-3 

A 2-1 

A 2-2 

A 2-3 

A 2-4 

B1-1 

20-30 

40-50 

49-59 

140-150 

160-169 

168-178 

181-190 

1-5 
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Table 2 

Atomic charges of amino acid residue His-41 of SARS CoV M pro , in 
atomic unit e 


Atom 

AMI 

Mull 

C J 

AMI 

q ESP 

HF/6-31G 

Mull 

HF/6-31G 

q ESP 

N 

-0.5314 

0.5569 

-0.9054 

-0.4429 

H(N) 

0.2987 

-0.5387 

0.4113 

0.1614 

C. 

-0.05334 

2.0820 

-0.0462 

0.1192 

H a 

0.2611 

-0.7777 

0.3081 

0.2461 

C 

0.3174 

-1.9381 

0.7638 

0.3254 

O 

-0.3257 

0.0303 

-0.5750 

-0.5518 

Cp 

-0.2302 

-1.2490 

-0.3088 

-0.1595 

Hp 3 

0.2633 

-0.1937 

0.3268 

0.2057 

H p2 

0.1628 

0.1624 

0.1903 

0.0576 

c y 

-0.0590 

1.7519 

0.2424 

0.2770 

N S1 

-0.2861 

-1.7978 

-0.8863 

-0.4094 

H 51 (+) 

0.4132 

0.7184 

0.4821 

0.4318 

Cs2 

-0.1669 

1.5835 

0.1266 

0.0852 

H 52 

0.3390 

-1.5772 

0.3350 

0.0992 

c el 

0.0476 

0.1000 

0.5776 

0.2349 

H s i 

0.3489 

1.6428 

0.3430 

0.2178 

n s2 

-0.2709 

-3.4887 

-0.9061 

-0.4587 

h s2 

0.4020 

2.1826 

0.46957 

0.4604 


from Table 2 that the negative charge of N§i (—0.4094) 
and the positive charge of (0.4318) are large in mag¬ 
nitude and are close to the peptide bond Gln-Ser. 
Therefore, the proton attracts the n electron density 
from the peptide bond Gln-Ser so as to weaken this 
chemical bond. The polar nitrogen N§i attracts one 


A 



hydrogen of bridge water and makes the water mol¬ 
ecule dissociated, causing the nucleophilic attack of hy¬ 
droxyl group OH - to the C(CO) and the electrophilic 
attack of H + to N(NH) on the two sides of peptide bond 
Gln-Ser. The sulfur atom S in Cys-145 forms a hydro¬ 
gen bond with H(NH) and helps the proteolytic reaction 
(Fig. 3B). 

Chemical modification of cleavable peptides 

The peptide bond is considered a pseudo n bond. The 
six atoms in the peptide bond (Gln)Ca CO-NHCa (Ser) 
are almost in a plane. Fig. 5A shows the electron density 
contour map of the peptide bond Gln-Ser in the gaseous 
phase on the 7i plane containing the carbonyl group CO 
of Gin and the NH group of Ser. The electron density 
contours surrounding atoms N(NH) and C(CO) form 
two triangles reminiscent of the hybrid orbits of the 
sp“ type, and the 2p z orbits, which are perpendicular 
to the plane and form a n\ bond system. In Table 3, 
we list the calculated atomic charges of the six atoms 
on both sides of the peptide bond Gln-Ser obtained 
from the ab initio HF/6-31G* calculations in the gas¬ 
eous phase with the background charges of SARS- 
CoV M pi °. The negative charge of N(NH) of serine in 
the octapeptide increases to —0.8612 in the background 
of proteinase charges from the original value of —0.8534 
in the gaseous phase. Similarly, the corresponding 


B 



Fig. 5. (A) Electron density contour map of the peptide bond Gln-Ser in the n plane containing the carbonyl C and O atoms of Gin and N(NH) of 
Ser, as calculated in the gaseous phase. (B) Contour map of the differential electron density of peptide bond Gln-Ser in octapeptide AVLQSGFR, 
obtained by subtracting the electron density in the gaseous phase from the electron density in the proteinase background and solvent water. 
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Table 3 


Atomic charges of relevant atoms in peptide bond Gln-Ser of octapeptide, in atomic unit e 


Phase 

Gin side 



Ser side 




Cot 

C(CO) 

O(CO) 

Cot 

N(NH) 

H(NH) 

g ESP (Gas) 

-0.0133 

0.8534 

-0.6252 

0.4124 

-0.8050 

0.3528 

? ESP (M pro ) 

0.0020 

0.8612 

-0.6219 

0.3842 

-0.8190 

0.3497 


positive charge of the carbonyl carbon C(CO) in gluta¬ 
mine of the octapeptide increases to 0.8190 from the 
original 0.8050. This is favorable to the electrophilic 
attack of OH - to C(CO) and the nucleophilic attack 
of H + to N(NH) in the proteolytic reaction. 

Fig. 5B is the contour map of differential electronic 
density of the peptide bond Gln-Ser in the octapeptide 
AVLQSGFR, obtained by subtracting the electron den¬ 
sity in the gaseous phase from the electron density in the 
background [23] of SARS CoV M pro and solvent water 
molecules. In Fig. 5B the bold gray line is the 0 value 
line where the electron densities in the gaseous phase 
and in the protease background are equal. The solid thin 
lines show the regions where the computed electron den¬ 
sities are greater in the SARS CoV M pro background 
than in the gaseous phase, and the dashed thin lines 
show the areas where computed electron densities are 
smaller in the proteinase background than in the gas¬ 
eous phase. We find that along the peptide bond be¬ 
tween (Gln)C-N(Ser) the electron densities increase on 
the N(Ser) side and decrease on the C(Gln) side. This 
change is favorable for a nucleophilic attack of anion 
OH - to (Gln)C and for an electrophilic attack of cation 
H + to N(Ser). 

SARS CoV M pro has very high cleavage specificity 
and is highly conserved in the coronavirus family [7]. 
The cleavage specificity of SARS CoV M pro concen¬ 
trates on the sites R2, Rl, and Rl'. The development 
of peptides as clinically useful drugs is limited by their 
poor stability and bioavailability, which is due in part 
to their inability to readily cross membrane barriers such 
as the intestinal and blood-brain barriers. Systematic 
chemical modification strategies that convert peptides 


into drugs are an attractive research topic in current 
medicinal chemistry [28]. For the peptide inhibitor of 
proteinase, the chemical modification to cleavable octa¬ 
peptide should focus on the scissile peptide bond be¬ 
tween Ri and Rg to be cleaved by SARS CoV M pro . 
A successful example shows that after the peptide bond 
CO=NH is replaced by reduced single bond CH 2 —NH, 
the modified octapeptide rennin inhibitors have shown 
greatly increased affinity for rennin and strong resis¬ 
tance to enzyme hydrolysis [29]. 

In this study we suggest a similar modification to 
octapeptide AVLQSGFR. In the peptide bond (Gin) 
CO=NH(Ser) the carboxyl group CO of glutamine on 
subsite Ri is replaced by CH? or CF 2 or the NH of ser¬ 
ine on subsite Rg is replaced by CH 2 ; then the n bond 
system is broken and the modified octapeptide AVLQ 
SGFR may become a stable inhibitor for SARS-CoV 
M pro and an effective drug candidate against SARS. 
Here we show the possibility of the chemical modifica¬ 
tion through computational modeling. 

Table 4 lists the atomic charges of the six atoms on 
both sides of the hybrid peptide bond Gln-Ser after 
the chemical modification. Comparing with Table 3, 
we find that after CO is replaced by CH 2 , the charge 
of the carbon atom in CH 2 turns negative (—0.6572) 
from originally a positive value (0.8612) in the CO 
group; hence the nucleophilic attack by OH - becomes 
impossible. Furthermore, the negative charge of N(NH) 
on the Ser side decreases to —0.2892 from —0.8190 in 
Table 3, and hence the electrophilic attack by H + also 
becomes more difficult. Likewise, the modifications 
CO -a CF 2 and NH -a CH 2 will make the proteolytic 
reaction more difficult also. 


Table 4 

Atomic charges in peptide bond Gln-Ser after modification, in atomic unit e 


Gin side Ser side 


CO —> ch 2 

Cot 

C(CH 2 ) 

H(CH 2 ) 

Cot 

N(NH) 

H(NH) 

^SF(Gas) 

0.0770 

-0.6756 

0.1351 

-0.1689 

-0.2473 

0.1910 

g ESP (M pr 0) 

0.0867 

-0.6572 

0.1251 

-0.2121 

-0.2892 

0.1807 

CO -a cf 2 

Cot 

C(CF 2 ) 

F(CF 2 ) 

Cot 

N(NH) 

H(NH) 

^SF(Gas) 

0.0814 

0.3679 

-0.0625 

0.2955 

-0.6854 

0.3143 

g ESP (M pr 0) 

0.1182 

0.3092 

-0.0653 

0.2597 

-0.6483 

0.3085 

NH -a CH 2 

Cot 

C(CO) 

O(CO) 

Cot 

C(CH 2 ) 

H(CH 2 ) 

^SF(Gas) 

0.0311 

0.8581 

-0.5980 

0.4387 

-0.6707 

0.1546 

g ESP (M pr 0) 

0.0551 

0.8437 

-0.5940 

0.4086 

-0.6460 

0.1511 
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Discussion 

Searching for cleavable peptides is the first step in 
finding potential inhibitors against the proteinase. Even 
within a limited range of octapeptides, it is by no means 
an easy job because the possible different patterns of 
octapeptides are as many as 20 8 = 2.56x 10 10 . However, 
some algorithms, such as the vectorized sequence-cou¬ 
pled algorithm [14] and the discriminant function algo¬ 
rithm [15], can be successfully used to narrow the 
search scope. Meanwhile, similar methods that can be 
used to search the cleavable peptides by SARS enzyme 
have been developed [5,6]. 

The octapeptide NH 2 -Ala-Val-Leu-Gln-Ser-Gly- 
Phe-Arg-COOH is the first designed on the basis of 
the molecular structure of SARS-CoV M pro [9] and it 
has been experimentally demonstrated to be a cleavable 
peptide with high bioactivity [12]. A detailed study of 
the cleavage mechanism needs a great deal of experimen¬ 
tal work. However, computational simulation may pro¬ 
vide useful insights for an in-depth understanding of the 
mechanism. The cleavable octapeptide could be changed 
into an effective inhibitor of SARS-CoV M pi ° or a drug 
candidate against SARS after a proper chemical modifi¬ 
cation to replace the cleavable peptide bond with a non- 
cleavable bond. The octapeptide thus modified will lose 
its susceptibility although it can still tightly bind to the 
active site and become a competent inhibitor, just like 
the role of a distorted key in blocking the access of other 
keys to a lock, as elucidated by Chou [15] in a review 
paper. 
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